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Abstract 

We discuss the fully non-linear formulation of multigravity. The concept of universality 
classes of effective Lagrangians describing bigravity, which is the simplest form of multi- 
gravity, is introduced. We show that non-linear multigravity theories can naturally arise in 
several different physical contexts: brane configurations, certain Kaluza-Klein reductions 
and some non-commutative geometry models. The formal and phenomenological aspects of 
multigravity (including the problems linked to the linearized theory of massive gravitons) 
are briefly discussed. 
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1 Introduction 



One of the most important problems which is facing theoretical physics now is the blending 
of the Standard Model (SM) with General Relativity (GR). Whatever way we choose (the 
most popular ones nowadays are based on some multidimensional constructions involving 
extended objects), nobody doubts that it will definitely modify physics at short scales. On 
the other hand, the current general paradigm is to keep General Relativity unchanged at 
large scales, but to add new forms of gravitating matter beyond the Standard Model (dark 
matter, dark energy) for explaining pressing astrophysical and cosmological facts such as 
galactic rotational curves and the accelerating universe. In the present paper, we consider 
an alternative paradigm: a modification of General Relativity at large scales as a possible 
explanation of some pressing cosmological issues (notably cosmic acceleration). 

The modification of GR that we are going to consider is linked to the issue of "massive 
gravity" (for very light gravitons, with Compton wavelength of cosmological scale). A 
generic prediction of multidimensional constructions is the existence of massive gravitons. 
In particular, any Kaluza-Klein (KK) model predicts, besides a massless graviton, the 
presence of an infinite tower of massive gravitons. However, it seems impossible to use the 
tower of massive KK gravitons to modify gravity at large scales. Indeed, its spectrum is 
generically regularly spaced (as illustrated on Fig. la), so that, even if the first mode were 
very light (i.e. of cosmological Compton wavelength), there would exist no regime where the 
first mode (or first few modes) would be important, and where one could truncate away the 
rest of the tower of massive states. In other words, as soon as the first mode is important, 
we open the extra KK dimensions (see, however, below). The situation is, however, different 
in some brane models. In particular, Refs. [l]-[4] discovered the possibility (illustrated in 
Fig. lb or Fig. lc) of having a hierarchical gap, mi <C vrt2, between the first mode (or 
first group, or even band, of modes) and the tower of higher modes. This situation, called 
multigravity (see [5] for a review and [6] for detailed presentation), makes it possible to 
envisage an effective four-dimensional theory which contains only the massless and ultra- 
light gravitons and discards the states of mass m > m 2 . The constructions [l]-[4] predict 
see-saw-like spectra, miMp^ ~ ^2 +7 > with 7 interpolating [5] between [1] and 1 [2]. 
Such spectra are naturally compatible with the phenomenologically interesting situation 
where m{ 1 is of cosmological order, while m^ 1 is smaller than the millimetre scale. 

So far multigravity was only analyzed in the linearized approximation. The main em- 
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a) b) c) 

Figure 1: Regular spectrum on (Fig.l a) versus bigravity (Fig.l b) or quazi- localized gravity 
(Fig. 1 c). The last spectrum is continuous but the first band is very narrow in comparison 
with the gap between bands. 

phasis of this paper is to provide a fully non-linear formulation of multigravity, i.e. to write 
down, and analyze, a class of consistent effective four-dimensional Lagrangians, describ- 
ing, in some limit, the light-mode truncation of the hierarchical spectra of Figs, lb or lc. 
Though we shall illustrate below our approach in the context of particular multidimensional 
realizations (notably brane models exhibiting multilocalization [1],[7] or quasi-localization 
[2], [8], [9]), we view our considerations as concerning a very general phenomenon: the con- 
cept of Weakly Coupled Worlds (WCW). The concept of WCW is very simple: one assumes 
that there are several Universes (labelled by % = 1,...,N), each endowed with its own 
metric and set of matter fields which are coupled only through some mixing 

of their gravitational fields. We require that the theory describing the WCW be near a 
point of enhanced symmetry, in the sense that there exists a limit (say as some parameter 
A — > 0) where the theory contains iV diffeomorphism-like symmetries, corresponding to iV 
massless gravitons. A recent theorem [10] has proven that the only consistent non-linear 
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theory involving N massless gravitons is the sum of N decoupled GR-type actions 

N 

So = £%*,<£>*], (i) 
i=i 

with (we use the signature — h ++) 

S[ 9i , $J = | rf 4 x ^g~ [M^R(gi) -A l + L( 9i , $,)] • (2) 

Therefore, the only consistent action for a theory of worlds coupled only through gravity is 
of the form 

N 

Stat = £ S[g i: $i] + A S mt (g 1 ,g 2 , . . . , g N ) ■ (3) 

i=i 

When A — > 0, the N worlds are non interacting (which implies that, from the point of view 
of any observer in one world, the other worlds have only a meta-physical existence), and the 
theory has the enormous symmetry rLDiff(j), where each diffeomorphism group Diff(j) acts 
separately on its own metric g^^ v and matter fields In the interacting case, A ^ 0, 

the symmetry of the full action must (again because of the theorem [10]) be reduced to 
(at most) one group of diffeomorphisms: the diagonal group of common diffeomorphisms 
transforming all metrics as 

5 g® = e X d x g® + d, e x g® + d u e x g® = D® e u + D® e, . (4) 

This symmetry restricts the interaction term A Sint^i, • • • , <7iv) to depend only on the invari- 
ants one can make with several metrics. This even leaves room for extra kinetic terms built 
from covariant derivatives such as g?X gffl (such terms do not exist in the case of one 

metric because D® g® = 0). However, in view of the many potential diseases associated to 
modifications of the standard Einsteinian kinetic terms, and in the spirit of describing the 
class of interaction terms most relevant at large scales, 3 i.e. containing the lowest possible 
number of derivatives ( namely zero, as expected from a generalization of the mass terms 
that appear in linearized multigravity), we shall only consider ultra-local interaction terms, 
i.e. 

AS , int = -/i 4 J d 4 xV{g!(x), . . . ,g N {x)) , (5) 

3 See Section 2.1 below for further discussion about extra kinetic terms. 
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where /i is a mass scale (henceforth replacing A as "small parameter") and where V is a 
scalar density made out of the values of the N metrics at the same "point" . We assume, 
for simplicity, that the N weakly coupled worlds "live" on the same abstract manifold, 
i.e., in other terms, that one is given a family of (smooth) canonical one-to-one maps: 
world(j) — > world(j). 

The aim of this paper is threefold: (i) to motivate the possibility of the effective action 
(3), (5) by considering several different specific models (brane models, Kaluza-Klein mod- 
els and non-commutative geometry ideas); (ii) to delineate and parametrize the various 
"universality classes" of non-linear multigravity; and (iii) to sketch the main qualitative 
consequences of such non-linear multigravity theories and to contrast them with the usual 
paradigm of "massless plus massive gravitons" which is based on a linearized approximation. 

It should be noted that theories defined by (3), (5) (in the "bigravity" case: N = 2) were 
first introduced in the seventies [11] as a model for describing a sector of hadronic physics 
where a massive spin-2 field (the "/ meson", with "Planck mass" Mf ~ 1 GeV in Eq. (2)) 
plays a dominant role. It was then called "strong gravity" or the "f-g theory" . Our work not 
only proposes to revive, within a new (purely "gravitational") physical context, this early 
proposal, but initiates the task of systematically studying the general phenomenological 
consequences of the action (3), (5). The present paper will only briefly sketch the new 
physical paradigm following from such actions. In subsequent papers, we shall discuss 
in detail the cosmological consequences of such theories [12], as well as its strong-field 
phenomenology [13]. 

2 Universality Classes of Bigravity Effective 
Lagrangians 

For simplicity, we focus, in this paper, on the case of "bigravity" , i.e. N = 2. Understanding 
this case is a prerequisite for understanding the general multigravity case (N > 2). Let us 
note also that the bigravity "potentials" that we discuss here can be immediately used in 
the general case. Indeed a rather general class of "iV-metric potentials" V(<?(i), . . . , <7(jv)), 
Eq. (5), is the class containing only "two-metric interactions": V = V(9(i), 9{j))- F° r 
instance, one can define a "crystal- like" many-world with "nearest neighbour" interactions 
only V = Z)j V(<7(j), <7(j+i)). It is interesting to note that the continuum limit (N — > oo) for 
some suitable "nearest neighbour" interactions can mimic the propogation of gravity in a 
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higher-dimensional space, i.e. the term / dya A (y) \J — g(x, y) [tr((/ 1 d y g) 2 — (tr g 1 d y g) 2 ] . 
as in (29) below. See [12] for further discussion of this subject. 

2.1 Parametrization of invariants 

Using, when N = 2, the notation g^ = g L (for "Left") and gp) = ga (for "Right"), and 

factoring a conventional "average volume factor" (giga) 1 ^ 4 out of the scalar density V 4 , 
the generic bigravity action reads 



Note that the bigravity action (6) contains 5 dimensionfull parameters: two "Planck masses" 
M L and M R , two cosmological constants A L , A R (with dimensions M 4 ), and the "coupling 
mass scale" fi. 

Before proceeding, we note that the mass scale fx, entering Eq. (6), will be treated here 
as a constant parameter determining the coupling of the two worlds. However, one should 
keep in mind the possibility that it be replaced by a fluctuating field. This is suggested, 
in particular, by the brane realizations of bigravity where the value of \i depends on the 
physical distance between the branes, which is controlled by dilaton/radion fields. A more 
general model where \i — > fJ>(x), and where one adds a kinetic term for fx(x), may play an 
important role in addressing crucial cosmological issues (such as inflation) in the context 
of multigravity theories. 

Before we shall proceed further let us make two additional comments 

• We shall treat V as a potential here, i.e. as an ultra-local function of gi and 
g R . As already mentioned above one could also include extra kinetic terms like 
gt iu D L,R g <Tf> LD L,R gL R ^ eta For examp i e m i X ed terms like y/-g L R(9R) or ^/-g R R(gL) 

are of these type. However, let us emphasise again that the fundamental concept of 
WCW explored here is that one is required to be near the point of enhanced symmetry 
4 We could, instead, have factored out of V the other natural symmetric density: y/—gL + \J—9r- 



S = j d A x^g- L (M 2 L R{g L )-A L )+ J ' d A x ^g~ L L($ L , g L ) + 



/ d A x^g^{M 2 R R{g R )-A R ) + J d 4 Xv ^L^ R ,g R ) 
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and according to the theorem proven in [10] in this point mixing of different metrics 
is forbidden. Because of this all higher derivative terms mixing different metrics must 
be supressed by the small parameter /z 4 and enter as // / (D 2 g / M 2 ) . One can see that 
in the long-wave limit D 2 g/M 2 ~ k 2 /M 2 « 1 one can neglect their contribution in 
comparison with a potential term. 

• Let us stress that we did not introduce any direct coupling between the two worlds 
exept the indirect one mediated by gravity. It means that all matter which can be 
observed by any observer is locally coupled only to one world or the other one. Thus 
any local experimental check of the equivalence principle is the same as in the General 
Relativity - all locally mutually observable matter moves in the same metric. One 
can consider other cases but such a study will be left for future publications. 

The common diffeomorphism invariance (4) restricts the scalar potential V{gL,gR) en- 
tering Eq. (6) to depend only on the invariants of the mixed tensor H = gj} g^, i.e. 

Hi = g^ gBav . (7) 

In 4 dimensions, there are (because of Cayley's theorem) only 4 independent scalar invari- 
ants which can be made from H. For instance, using a matrix notation for H, one can take 
the first 4 traces of the powers of the matrix H, say 

r n = tr(H n ); n=l,2,3,4. (8) 

Let us introduce the 4 eigenvalues A a (a = 0, . . . , 3) of H£, i.e. the 4 eigenvalues of the 
metric g^ with respect to g^. They can be defined either by r n = S a A", or by writing the 
two metrics in a special bi-orthogonal vierbein such that 

o L = - e ° e° + e 1 e 1 + e 2 e 2 + e 3 e 3 

9% = -Ao e° el + A x e\ e\ + A 2 e 2 e 2 + A 3 ej e 3 . (9) 

It is easily seen that, apart from an exceptional case (where two eigenvalues coincide, 
and correspond to a null eigenvector), it is generically possible to write Eq. (9), though 
maybe with a complex vierbein e£. Indeed, two (but at most two) eigenvalues, say Ao, Ai 
(one of which necessarily corresponds to a time-like direction) can become complex. We 
shall focus on the case where the 4 eigenvalues A , Ai, A 2 , A 3 are real and positive. As 
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we shall only deal with symmetric functions of the eigenvalues, this restriction is mainly a 
notational convenience which can be relaxed by analytic continuation. It is then convenient 
to parametrize the invariants of H = gj} g R by means of the logarithms of the eigenvalues 
of if: 

H a = \n\ a - X a = e" a , (10) 

(the /i a 's should not be confused with the mass scale /z 4 in front of V) and to introduce, as 
basis of independent scalars, the 4 symmetric polynomials 

^ = £/C (ii) 

a 

With this notation, our first result is that the most general (densitized) potential can 
be written as 

^y{9L,9R) = ^ 4 {9l9r)^ V(a 1 ,a 2 ,a 3 ,a 4 ) , (12) 
where V is an arbitrary function of the 4 cr a 's. 

2.2 Universality classes 

In the same way as the various mathematical forms of the Landau free energy define uni- 
versality classes of phase transitions, we can define universality classes of bigravity theories 
by considering as equivalent the functions V(a a ) leading to (essentially) the same multi- 
gravitational phenomenology. As we shall see below and in [12, 13], some of the important 
qualitative features of the function V(a a ) are: (i) its behaviour near a a — > 0, (ii) its be- 
haviour when a a — > oo, and (iii) the existence or non-existence of "critical points" where 
some derivatives dV/da a vanish. 

As a first example of a universality class, we can define the class of Vs which reduce (in 
absence of cosmological constants in (6), i.e. Aj, = A^ = 0), in the linearized approximation, 
to the Pauli-Fierz mass term ~ h^h? — hft h v v . The linearized approximation corresponds 
to the particular case where gi,^ arid g Rflu are both near the same flat metric rj^, i.e. 
9l^ = r\ixv + h% v , g Rflu = rj^ + h% v , with h L < 1 and h R < 1. In this limit the above object 

H = g^ 1 g R reads Hjf ~ 5% + h Rl/ — h? Lv (where the indices on h R and hL are raised by rf v ). 
It is then seen that the eigenvalues of H are A a ~ 1 + /i a , where /i fl < 1 are the eigenvalues 
of h Rv — h1 u . With the identification of the massive graviton mode h^ as h^ v = h^ u — h^ v 
(see below) , one then sees that the Pauli-Fierz mass term ~ h£h' ' x — h^ h v v is obtained if the 



7 



function V(a a ) behaves (modulo a positive factor that can be absorbed in the mass scale 
/i 4 ) as 

V(<j a ) ~ (7 2 - of = YsVl ~ ( 2 AM , when/i a ^0. (13) 

a \ a / 

The behaviour (13) near /x a — ► defines the universality class of "Pauli-Fierz-like" bigravity. 
Note that one can imagine a case where the potential V does not have quadratic terms 
when ii a — > 0. In the linearized approximation, one would see two massless gravitons, while 
the full theory would contain two interacting metric field gi,gR (and only one common 
diffeomorphism invariance). 

As a second example of the concept of universality class, we can define the class of 
potentials V(gL,gii) which are symmetric under the exchange g^ <-> g^. It is easily seen 
that under the exchange L <-> R, the eigenvalues A a get inverted (A a — > A" 1 ) so that 
the logarithmic eigenvalues /i a change sign: fj, a — > —fJ> a - The class of exchange-symmetric 
potentials therefore corresponds to the class of functions V(/i a ) which are even in the // a 's. 
In terms of the cr's this becomes V(o~i, <7 2 , cr 3 , <r 4 ) = V(— a±, <7 2 , — <7 3 , cr 4 ). 

As a further example of universality class, we can consider the class of functions which 
depend only on the first two invariants a\ = I] a /i a and cr 2 = J^a^l '■ V — V(°~i,°~2)- 
We shall see that this class appears naturally in brane models, and our (preliminary) 
investigations suggest that this class might be general enough to describe all the possible 
qualitative features of a general bigravity theory. 

2.3 Equations of motion 

The equations of motion derived from the bigravity action read 

2 Ml - 1 R(g L )) + Al g£, = t\ v + r£ , 

2 M\ (R»M R ) - \ 9* R(9 R )) + A fl g% = t% + T£ . (14) 

Here T^" = 2(— <7l) -1 / 2 5 S^ aUcr /S g^ u denotes the stress-energy tensor of the matter on the 

left brane (S^ atter = J d 4 x y/-g L L($ L , 9l)), while = gT 9l denotes the effective 
stress-energy tensor (as seen on the left brane) associated to the coupling term S- mt = 
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2 5 S- mt (g L , g R ) 



Sg L 



= -2/x 4 



9r 
v 9l 



,xv y . dV(g L ,g R ) 
9l J + 



(15) 



The corresponding expressions for the right brane are obtained by the exchange L <-> i?. 
For instance 



^ ^ , dV(g L ,g R ) 
9R 4 + 9< 



(16) 



The Bianchi identities [R^u ~ \R L d^u) = o)> an d the conservation of the material 
energy tensor [D V L T^ v = 0; when the matter equations of motion are satisfied) imply the 
constraints: 

0. (17) 



and D R t R v 



Actually these two constraints are not independent because the invariance of S[ nt under the 
unbroken diagonal diffeomorphism group implies the identity 

V^glDl t L ^ u + ^g~ R D R t% = . 

The explicit expressions of the derivative terms dV/dg^ R in Eqs. (15), (16) tends to be 
rather complicated. However, they acquire a simple form when written in the special 
frames with respect to which both and g R v are diagonalized (such as in Eq. (9)). The 
mixed components of vL, and t R v with respect to any such frame (which can differ from the 
particular of (9) by arbitrary rescalings — > ( a e" because such rescalings leave t a La and 
t Ra invariants) take the simple form: (no summation on the frame index a) 



t a 



-2/i 4 e^ 



V 



'V | dV" 
4 dfi a/ 



(18) 



with vanishing of the off-diagonal components (we recall: a\ = X^A^)- 

Here, we considered the scalar potential as a function of the /x a 's. If V is given as a 
function of the cr n 's, Eq. (11), the derivative entering Eqs. (18) takes the explicit: 



dV(ai, . . . ,a 4 ) 



dV dV 2 dV A ,0V 



da 2 



da 3 



da 4 



(19) 



This explicit expression illustrates the third type ((iii)) of universality class mentioned 
above: If there exist "critical points" where d a2 V — d a3 V — d a4 V — (without restriction 
on d ai V), such points give rise to a t^ Lu and a t Ru with the local "equation of state" 
^lo = = = ^L3 ( an d similarly for t R ), i.e. such that oc g^ u and oc gjf v . In some 
cases, such critical points can be "fixed points" and can give rise (in the "vacuum case", 
i.e. in absence of "material" T^ R ) to bi-(A)dS solutions of the coupled field equations. 
Note in this respect that the "perturbative limit" fi a = is a critical point in the sense 
that <9 Ma V = d ai V, independently of the value of a, so that /i a — (i.e. g^ v = g^ u ) can be 
a (perturbative) fixed point of the coupled vacuum equations, corresponding to a bi-(A)dS 
solution, if the corresponding (constant) curvature A (Rl v = = R Ru ) satisfies the two 
equations 

—2 A M| + A L = -2/i 4 



-2 A M R + A R = -2 fi 4 



V 



IJ.a=0 



Ma=0 



(20) 



In the "Pauli-Fierz" universality class the right-hand sides of Eqs. (20) vanish and one has 
the usual relation A = A L /(2M|) (with the constraint A L jM\ = A R /M R ). In more general 
classes the coupling between the two worlds can modify the usual link between A and Al,r- 



2.4 Single "massive graviton" as a limiting case of bigravity 

Let us consider the formal limit M R — > oo in the action (6) (and the field equations 
(14)). In this limit the metric g R is (formally) frozen into some given "background" metric 
G^u '■ 9r —> G, with G^ solution of R 11V {G) = XG^, where A = lim(A^/(2 M R )) can be 
zero, or can be arranged to take any fixed real value. This leaves us with an action for a 
single dynamical metric g = g^ of the form 

S = Jd i x^(MlR(g)-A L )-^Jd 4 x(gG^V(G- 1 g). (21) 

If V belongs to the Pauli-Fierz universality class, g^ u = G^ u is a solution of the equations of 
motion (if A = A L /(2 M|)), and the small excitations of g^ around G^ u describe a "massive 
graviton" (propagating in an Einstein space). But the behaviour of the large excitations of 
g^v are described by the non-linear action (21) instead of the usual quadratic Pauli-Fierz 
action. 
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The action (21) is (formally) generally covariant: when g is transformed as (4), the 
frozen metric G must also be transformed as 5G^ = D*? e v + -D^e M . These fluctuations 
of G (which do not change the background curvature invariants) are playing the same role 
as the Goldstone degrees of freedom in the Higgs mechanism for gauge fields. In a recent 
paper [14] these Goldstone degrees of freedom were discussed for the single AdS4 brane case, 
using an holographic description of five-dimensional gravity in terms of a four-dimensional 
CFT, and it was shown that there is indeed a vector field which provides extra components 
to the graviton. 

3 Specific Examples of Bigravity Effective 
Lagrangians 

After having discussed general possible structural features of bigravity effective Lagrangians, 
we shall consider specific physical models in which such Lagrangians arise. We consider in 
turn: (i) brane models, (ii) Kaluza-Klein models, and (iii) non-commutative-type models. 
Beforehand let us note that the work in the seventies that first considered bigravity models 
did not have any underlying physical models from which they could derive some specific 
potentials V(gL, gn). They made up some non-linear generalizations of the quadratic Pauli- 
Fierz mass term. For instance, they particularly considered the one-parameter family of 
models with 

V(g L ,g R ) oc {^j (<4«4 - g L a ,g^)(g a / - g?)(<% - 9?) . (22) 
3.1 Brane models 

Let us start by briefly recalling why (multi-)brane models naturally give rise to "multi- 
gravity". For more details the reader is advised to look at the original papers, and/or 
at reviews such as, [15], [5]. Before explaining how several worlds can be gravitationally 
"weakly coupled" , let us recall that the paradigmatic brane example of a separate (gravita- 
tionally decoupled) brane world is a Randall-Sundrum (RS) scenario, i.e. a flat 3-brane in 
AdS(5), with jump conditions on the brane (coming from an assumed Z2 symmetry) able 
to "localize" the 5-dimensional graviton as a massless excitation propagating (as a "surface 
wave") in the vicinity of the brane [16], [17]. Putting the brane at the point y = (where 
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y is the "fifth", transversal coordinate, and where one requires the Z 2 symmetry y — > —y), 
the background 5-dimensional geometry is (see Fig. 2a) 

ds 2 = e- 2a ^r]^dx' x dx v + dy 2 = e~ 2a{z) [ Vliu dx^dx u + dz 2 ] . (23) 

The warp factor behaves as a(y) ~ \y\ ~ ln(l + \z\). The fluctuations near the background 
metric are studied by writing: 

ds 2 = e~ 2cTiz) [(77^ + h^(x, z)) dx li dx v + dz 2 ] . (24) 

The field h^ u (x,z) is expanded in terms of the graviton and KK plane wave states : 
h^ix.z) = J2^=o exp (§ o") h^J(x)^ n (z) where the exp (| o^j factor in the expansion is 
necessary for the functions ^ n \z) to obey an ordinary Schrodinger equation: 

{-d 2 + V(z)}y n (z)=m 2 n V n (y). (25) 

Here the potential V(z) = [dAj dz) 2 — d 2 Aj dz 2 where A = 3a(z)/2. Qualitatively it is made 
up of an attractive ^-function potential plus a smoothing term (due to the AdS geometry) 
that gives the attractive potentials a "volcano" form. An interesting characteristic of this 
potential is that it gives rise to a (massless) normalizable zero mode 

= exp[-A(z)] = exp[~ a{z)\ (26) 

One can show (see for details [18] and references therein) that the normalization factor 
/ dz^Q = J dz exp[—3a(z)} also relates the fundamental five-dimensional mass scale M 5 to 
the four-dimensional Planck mass M p , namely M 2 = Mf / dzexp[— 3a(z)]. 

One can consider now multibrane configurations where the warped metric is a bounce 
as on Fig. 2b. By analysing the spectrum in this case one can easily see that in the case of 
an infinitely large separation between the branes massless gravitons are localized on both 
of them. But then, according to basic properties of the Shrodinger equation, when the 
separation is finite, the degeneracy between the two massless modes is removed and one 
ends up with one massless and one ultralight massive graviton. The prototype model of this 

class was the "H h bigravity" model [1] with two positive tension flat branes ("+" branes) 

separated at the bounce position L\ by one intermediate negative tension flat brane ("— " 
brane) in an AdS$ bulk. The task of finding the KK spectrum reduces to a simple quantum 
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e -2G(y) 



,- 2o(y) 




a) 




b) 



Figure 2: Warped metric for single flat brane (Fig. 2a) and bounce for a two-brane configu- 
ration (Fig. 2b). L 2 is the separation between branes and L\ is the position of the bounce. 
If Li <C L 2 the metric is mostly concentrated on a right brane and if L 2 — L\ <C Li then it 
is concentrated on the left brane. 



mechanical problem. It is simple to see that the model (as every compact model) has a 
massless graviton that corresponds to the ground state of the system whose wave function 
follows the "warp" factor. Then it is easy to see that (say, for simplicity, in the symmetric 
case L 2 = 2L ± ) there should be a state with wave function antisymmetric with respect to 
the minimum of the "warp" factor, whose mass splitting from the massless graviton will 
be very small compared to the masses of the higher levels. Because the "warp" factor is 
exponential the difference in mass behaviour of the first and the rest of the KK states is also 
exponential. This allowed for the construction of a linear "bigravity" model in which the 
remainder of the KK tower does not affect gravity beyond the millimetre bound. Soon after 
this model, other models were discussed. Some of them, for example the "quasi-localized" 

GRS model [2] and a more general H h multigravity model [3] , [4] also used dynamical 

negative tension branes. In other models, like the +- 1- model with two AdS 4 branes [19] (or 
the limiting case of one single AdS^ brane when the second one is moved to infinity [20]) 
or in a six-dimensional case [21], there are no negative tension branes and so no problems 
emerge with ghost-like radion states. Models with moving branes in which one also can 
get "warped" factors were discussed in [22]- [23]. Finally there is a whole zoo of different 
models in which one can get modification of gravity at large scales. 

To be specific, let us consider the H h model and let us now derive the fully non-linear 

bigravity action it gives rise to in the weak-coupling limit where we keep only the dominant 
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terms in the exponentially small ("tunnelling") coupling between the two positive tension 
branes. We are going to ignore the fact that there is a ghost-like radion field due to the 
existence of the negative tension brane in this model [24]. Anyway we freeze all dilaton and 
radion degrees of freedom. 

The action describing the full 5-dimensional configuration is (in units where the five- 
dimensional Planck mass is set to one) 



S = J d A x J dy 



R(G)-V b -J2V l 5(y-y t ) 



(27) 



Here, G denotes the 5-dimensional metric, Vb the bulk cosmo logical constant, and Vi the 
tensions of the branes (the index i takes, in our case, three values corresponding to the three 
branes: e.g. V± — Vl > 0, V 2 < and V3 = Vr > 0). We generalize the linear fluctuation 
ansatz (24) by writing the 5-dimensional metric as (/i, v — 0, 1, 2, 3) 

dsl = a 2 {y) g^{x, y) dx» dx v + dy 2 , (28) 

where a 2 (y) = e~ 2cr( ^ is the background warp factor. We assume here that the degrees 
of freedom associated to the fluctuations of: (i) the warp factor ("dilaton"), and (ii) the 
distance between the branes ("radions") are all frozen. The detailed mechanism of how 
to do that is not important for us now. For example, one can add extra terms in the 
action (27) that give large enough mass terms to these fluctuations (say with submillimeter 
Compton wave length) following [25] (of course for those radion fields which are ghost-like 
one has to add tachyonic mass terms). The fluctuations of the mixed components g m can be 
consistently set to zero, because of the Z 2 symmetry requirement. Inserting the ansatz (28) 
into the action (27) yields, after integration by parts and use of the background equations 
of motion for the warp factor (which allow one to dispose of all terms containing y-gradients 
of a(y)) the following action for 57^(2:, y) 

S = Jd'xjdy ^-g{x,y)[a 2 {y) R(g(x, y)) - I a\y) [tr(<T 1 d y g) 2 - (tr g' 1 d y gf] } . 

(29) 

Here, tr h 2 — (tr h) 2 = h% h 1 ^ — (hfy 2 where h% = (g^ 1 d y g)% = g^° d y g av . Note that this exact 
(after freezing the dilaton and the radions) action for the nonlinear dynamics of g^x^y) 
is still 5-dimensional. Note also that all explicit coupling to the branes have disappeared 
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(thanks to integration by parts). The crucial feature of (29) for our discussion of an effective 
4-dimensional bigravity action is the presence of the warp-factor dependent coefficients 
a 2 (y) and a A (y). It is the fact that these factors are exponentially localized on the two 
positive-tension branes (as shown in Fig. 2b which plots a 2 (y) = e~ 2a ^) which will allow 
for the derivation of an approximate 4-dimensional action. Though Eq. (29) was derived 
from an explicit 5-dimensional model, we expect that the general structure of Eq. (29), 
namely to have a curvature term with a weight function (here a 2 (y) dy) which localizes it 
on some branes, and a transverse-gradient term which also comes with a similarly localized 
weight function (here a A (y) dy), will hold in more general situations, like, for instance, the 
6-dimensional model of [21] (which is free from negative tension branes). Probably, in the 
latter model, if we assume that the excitations related to gradients in the sixth "angular 
type" direction 9 are frozen (i.e. massive enough), we shall get an effective action of the 
type (29) but, possibly, with weight factors which are somewhat modified (by the y- varying 
volume of the sixth circular dimension). To enhance the generality of our discussion, and 
cover such cases, we shall henceforth work with an action of the form (29) but with the 
replacements a 2 {y) — > a 2 (y), a A {y) — > a 4 (y), where a 2 {y) and a A {y) are two (unrelated) 
"weight" functions which are strongly localized around two branes. The generalization to 
the case of N branes is obvious. The essential features of a 2 (y) and a 4 (y) that will be needed 
in the following is that they are both positive and that: (i) a 2 (y) reaches maxima which 
are sharply localized on two branes, while (ii) a~l l {y) reaches a sharp maximum somewhere 
between the two branes. The crucial point is to realize that these generic conditions imply 
the following specific y-dependence of g^ v {x,y): as a function of y, g(y) is nearly constant 
everywhere, except in a "transition layer", located around the minimum of a±{y), where 
g(y) has a fast variation with y. In other words g(y) is a smoothed version of a Heaviside 
step function: g(y) ~ g 1 ^(y* — y)+g2 6{y — V*) where is the location of the minimum of 
di{y) and where g\ and g 2 are two different asymptotic values (which depend on x M when 
putting back everywhere the x-dependence). It is this transition-layer behaviour which 
allows us to derive an approximate 4-dimensional action for gi(x), g2{x). 

To understand intuitively this transition-layer behaviour we can assume that we normal- 
ize a^ 1 (y) so that it takes the value aj 1 ^*) = 1 at its maximum and then decreases to very 
small values as y gets away from y* (either way). Let us view the action (29) (with a 2 — > a 2 , 
a 4 — > a 4 ) as a "mechanical" Lagrangian for the motion of the particle g, when thinking of 
y as being "time". The "kinetic" terms are the last two terms quadratic in y-derivatives. 
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We then view a±(y) as the "mass" of the g-particle. This mass is of order unity around 
y = y*, and then increases to very large values on both sides. In other words, the g-particle 
is extremely "heavy" everywhere away from y*, and becomes relatively "light" only around 
y*, which makes it clear that g(y) will "move" very little away from y*, and that all "y- 
motion" will take place only around y*. Another way of seeing that what is important is 
to have separate maxima in a 2 (y) and a^ 1 (y) would be to consider the y-Hamiltonian, (in 
terms of the y-momentum rr = dC/d(d y g)) which is of the symbolic form: a 2 R + aj 1 n 2 . 
One can technically analyze the behaviour of g(y) in the transition-layer by zooming on the 
exact solution of the only relevant part of the "dynamics" near y = y*, namely the "kinetic 
terms" ~ a^(d y g) 2 ~ aj 1 it 2 . Note that a 2 takes very small values around y*, so that we 

can, in first approximation, neglect a 2 R with respect to aj 1 n 2 . This can be done exactly 
by changing the "time variable" . Indeed, in terms of the new "time" t, defined by 

dt = a^(y)dy, (30) 

the "kinetic" part of the action (29) reads simply 

S k = -Jd 4 xl, (31) 

where 



det g(t) ^{g^gy-itTg^gf . (32) 



Here g = dg/dt, and we leave implicit the x-dependence of g. Actually, (31) does not couple 
anymore t- and x-derivatives. Therefore we can solve the equations of motion derived from 
(31) separately for each point x, i.e. it is enough to solve (32) at each x. The action (32) 
is still a very non-linear action for the t-dynamics of a 4 x 4 matrix g^ v {t). However, it 
is exactly integrable. This is seen by exploiting the symmetries of /: (i) invariance under 
rigid SL(4) transformations of g^, and (ii) invariance under time translations. Note that 
we only have an SL(4) symmetry, and not a GL(4) one because of the presence of det g. In 
other words, the action is invariant under g^ v — > g' = g a p only when det A = 1. The 

first symmetry leads to the traceless mixed tensor constant of motion (in any spacetime 
dimension D) 

c^<-^KK- (33) 
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Here 7r^ = g va it^ where ir^ u is the "momentum" conjugate to g^ u , namely 



w = 



SI 



= y/=g{K^-KgT) 



(34) 



where iC" = g^ a g u/3 K af3 with K a p = | g a p being the usual "second fundamental form" . 
The second symmetry leads to the constancy of the "energy" : 



E = \ V=g~ [trig' 1 gf - (tr g' 1 gf] = ^g {K K; - K 2 ) 



(35) 



Contrary to what happens in the well-known Kasner solutions, we are not restricted here to 
the "zero-energy" shell (because of the influence of the curvature term a 2 R which changes 
the asymptotic behaviour of g on both sides of the "transition layer" that we are currently 
zooming into). This implies that the exact solution g^ v (t) is different from, and more 
complicated than, a Kasner solution. 

The exact solution is obtained by decomposing g^ v in its determinant (or better w = 
\J— det g) and its unimodular part, say 7^/ (— det g) 1 ^ '. Eq. (33) simply says that w 7 _1 7 
is the constant matrix 2 Off. This is immediately integrated to the matrix equation 



7 (t) = 7 (t )exp f2Cjf 



* dt 



w(t) i 



(36) 



To complete the solution for g(t) we need to know how its determinant — det g = w 2 depends 
on t. This is obtained by combining (33) with (35). This yields a first order differential 
equation for w(t): 

D - 1 



D 



w 2 = c 2 -Ew; c 2 =tiC 2 = C^C^. 



(37) 



In terms of the new parameter 



1 

x = - 

c 



° —t — B 



D - 1 2 



(3* 



where B is a constant of integration we get the solution 



w 



= ^J-detg = ^(l-x 2 ), 
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dt 



2 Id -i dx 



(39) 



w 



c V D 1-x 



This allows one to express the matrix g as an explicit function of x: 




) 



2_ 

D 



D - 1 



Jx 




9(0) exp 4 



D 



(40) 



where the matrix C is C% = c 1 C/f, i.e. is normalized so that tr C 2 = 1. 



The above exact solution for g(x), i.e. g(t), using (38), does not seem to involve any 
transition-layer behaviour. The transition-layer behaviour appears when we express g in 
terms of the original transverse variable y (which is the proper distance orthogonally to the 
branes). Indeed, when (qualitatively) integrating Eq. (30) to express t as a function of y, the 
sharp maximum of a^ x (y) around means that t(y) behaves essentially as a (smoothed) 
step function t(y) ~ t\ 9(y* — y) + t 2 6(y — y*). Inserting this sharp-transition behaviour 
into the smooth solution (g(x(t)) (40) then leads to the announced (smoothed) step-like 
behaviour of g(y), with the bonus that we now have in hand the (rather complicated) 
precise manner in which g(y) sharply (but smoothly) evolves in the transition region. It 
is interesting to make the link between the nonlinear transition of g(x,y) between the two 
positive-tension branes (which is a smoothed version of g(x, y) ~ gi(x) 0(y*—y)+g2(x) 6{y — 
y*)) and the result of linearized fluctuations which, as recalled above, is expressed as 



One indeed finds, when looking at the explicit results for the various mode functions 
ty n (y) that the first two modes (n = 0, 1; corresponding to the massless mode, and the 



e(x) = sign(x). Keeping only the first two modes is then equivalent to considering metric 
fluctuations of the form g^x, y) ~ i]^ u + hffi(x) +e(y — y*) hfy(x), which is fully consistent 
with our result for the fully nonlinear metric g(x,y) interpolating between a gi(x) and a 
g2{x) through a transition layer. When going beyond the step-function approximation, one 
can also check that the nontrivial transition behaviour (40) does also correspond (when lin- 
earized in g^u — iliiv) to a zoom on the (large k£) limit of the first mode e^ a ^\{y) (considered 
as a smoothed version of e(y — y*)). Note also that the characteristic width of the transition 




lightest mode) behave 




where 
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layer is Ay ~ k 1 , where k is the usual bulk curvature parameter defined such that the 
background solution has (d y a) 2 = k 2 (outside the branes), so that a(y) = k \y — yi\ near 
brane 1 (and a 2 (y) = e~ 2a = e 2fc ' 27 2/1 ' ) . There is a clean separation between the transition 
layer (around which is the location of the middle negative-tension brane) and the "lo- 
calization layers" (around y\ and y 2 , i.e. the locations of the two positive-tension branes), 
when kl ^> 1, where I denotes the smallest interbrane distance: I = mm(Li, L 2 — Li); see 
Fig. 2b. Because of the exponential dependence of the warp factors (and therefore of a 2 (y) 

and 04(1/) in the H h model), even a moderately large value of k£ suffices to ensure that 

the above (nonlinear) transition-layer approximation is valid up to exponentially smaller 
corrections. 

The exact, nonlinear transition-layer solution (40) interpolates between a certain metric 
gi(x) = g(x,yi) on the first brane, and another one g 2 (x) = g(x,y 2 ) on the second brane. 
Instead of viewing the exact solution (40) as the solution of a Cauchy problem (e.g. for given 
g 1 and gi), we should reexpress it as the solution of a "Lagrange- Feynman" problem, i.e. as 
the unique extremizing solution of the action (31), (32), for given "initial" and "final" values 
of g{y): i.e. for given gi(x) = g(x,yi) and g 2 (x) = g(x,y 2 ). We can also think of (32) as 
defining a certain Riemannian metric in the space of metrics g^ u . We are then considering 
the "geodesic" connecting some given initial point g\ to some given final point g 2 . Let 
g gi , g2 (y) denote this unique (parametrized) geodesic. The analysis above then leads us to 
estimate that a good approximation (when k£ ^> 1) to the effective action describing the 
dynamics of gi(x) and g 2 (x) is obtained by inserting the "geodesic" g gug2 (x,y) (computed 
for each point x) in the original full action (29), so that 



where lZ[g] = \J—g R(g), and where I[g\,g 2 ] is the value of the "geodesic" action (32) 
evaluated for the extremizing solution g gi , g2 (y) and integrated between yi and y 2 . The 
factors 2 in (42) come from the fact that we are assuming periodicity over y varying between 
yi and 2y 2 — y\- Calculating I[gi,g 2 ] from the exact solution (40) is somewhat complicated. 




(41) 



where (suppressing the x-dependence to focus on the y- dependence) 



rV2 

£[91,92] = 2 / dya 2 (y)K[g gug2 (y)] -2i"[#i,# 2 ] , 



(42) 
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Let us only give the final result (which is simpler than the necessary intermediate steps): 



where 



with 



I\9i,92] = (gig2)*V b (g 1 ,g 2 ), (43) 



V b (gi,g 2 ) = ^ (cosh/3 - cosh a) (44) 



py2 ry2 

T= dt= al\y)dy, (45) 



^W^TV^-^^V^T^' (46) 

<x=\<ti. (47) 

As above, o\ = T, a fi a , a 2 = S a/ u^ where \i a denote the logarithms of the eigenvalues of the 
matrix gf 1 g 2 , i.e. <j\ = tr In g^ 1 g 2 and o 2 = tr (In gf 1 g 2 ) 2 . The combination a 2 is £ a /I^ 
where /2 a = fi a — a\/D denote the logarithms of the eigenvalues of the unimodular metric 
ji 1 72, i.e. a 2 = tr (In 7-f 1 72) 2 . For added generality, we have left the dependence upon the 
brane (spacetime) dimension, though we have in mind here only D = 4. The weak-coupling 
parameter appearing in front of the interaction term / is the inverse of the total "t-time" 

T = J dt needed to interpolate between g\ and g 2 . We recall that, in the H V model, we 

have 04 = a 4 . An explicit computation then yields 

4ke 

T = w (48 > 

This exponentially large value (due to the exponentially small value of a 4 near the inter- 
mediate brane) corresponds to the expected exponentially small coupling between the two 
metrics on the positive-tension branes. 

To get an explicit bigravity action, one still needs to evaluate the first contribution in 
the Lagrangian (42). Neglecting exponentially small fractional contribution it is clear (in 
view of the localized behaviour of a 2 (y) and of the near y-constancy of g gi>g2 (y) outside of 
the transition-layer) that this contribution is well approximated by replacing g gi>g2 (y) by 
its (relevant) boundary value g\ or g 2 . Finally, the full brane-derived bigravity effective 
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Lagrangian density (in units where the coefficient of 1Z in the higher-dimensional theory is 
set to one) is 

£(<?!, g 2 ) = A 1 K( 9l ) + A 2 K(g 2 ) - 2I( 9l , g 2 ) , 

where 7Z(g) = y/^g~R(g), where (see Fig. 2b; we assume that y varies over a full period 
i-L 2 ,+L 2 \) 

/Li f2L,2—L\ 
dya 2 (y), A 2 = / dy a 2 (y) , (49) 

and where the "potential" term I(gi,g 2 ) is given by Eqs. (43), (44) above. 

It is easily checked that the potential (44) has the Pauli-Fierz limiting behaviour (13) in 
the limit fi a — > 0. One can then compute the corresponding Pauli-Fierz mass. One finds (in 
the symmetric case L 2 = 2L\ = 2£, for simplicity) m| F = 4T~ 1 A{ 1 . The explicit value, 

in the H h model, of A\ is A\ ~ 1/A;, so that we get mpF = 2\f2e~ 2kl k^ in agreement 

with the direct analysis of linearized fluctuations [1]. In the Appendix we further compare 
the nonlinear bigravity action to the linearized bigravity results already derived in the 
literature. In particular we check that they are fully consistent, even in the asymetric case 
L 2 ^2L V 

A full justification of the effective action (42) can, in principle, be obtained by explicitly 
considering the effect of corrections to our approximation g(x, y) ~ g gitg2 (x, y). For instance, 
we can write g{y) = g git g 2 {y) +C(y) where the correction £(y) vanishes, by definition, when 
y = y\ and y = y 2 . We can then expand £(y) = S n £ n sin 2nn r(y) where r = t/T varies 
(by definition) in the interval (0,1). [The condensed notation £ n denotes some C, n ^u( x )]- 
An analysis of the full action (expanded quadratically in the £'s), containing not only the 
"light fields" gi(x), g 2 (x), but the tower of "heavy fields" £„(x), shows that the mass of 
the heavy fields scale like ~ e~ M k, which is exponentially heavier (by a factor e +kt ) 
than the Pauli-Fierz mass scale. This confirms that the nonlinear bigravity action (42) is a 
good effective description when one considers configurations gi{x), g 2 (x) where the relevant 
gradients are small compared to m^. 

3.2 Kaluza-Klein Models 

As said in the Introduction, and sketched in Fig. 1, one expects generic Kaluza-Klein models 
to give rise to "regular spectra" containing no gap allowing one to separate a finite number 
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of light gravitons from an infinite tower of heavy ones. We wish, however, to emphasize the 
existence of a class of KK models where such a gap can exist. 

By KK model, we mean a higher- dimensional background geometry which decomposes 
as a direct (unwarped) product, ds 2 ot = ds 2 D + da 2 where ds 2 D = gffi(x) dx^ dx u (with, e.g., 
g$(x) = r]^ v ) and da 2 = jabiv) dy a dy b . When decomposing the fluctuations of the higher- 
dimensional metric g M N{x,y) = (g ixvi g^a,} gab) into representations of the symmetry group 
of gffl (say gffi = i]^ u ) one generally expects the squared mass spectrum of tensor (spin 
2) fluctuations 5 g^ u to be given by the spectrum of the scalar Laplacian on the (compact) 
internal manifold, say T, with metric da 2 = r y a b(y) dy a dy b . Let \ n > 0, with n — 0, 1, 2, . . . 
denote the latter spectrum, i.e. 7 -1 / 2 <9 a (7 1 / 2 / y ab d b ) 4> n (y) = —\ n (f) n (y). There is always 
a zero- mode, (f>o(y) = const., corresponding to A = 0, i.e. to a massless graviton. The 
question of the existence of a hierarchy allowing one to consider, for instance, an effective 
theory containing only the massless graviton and a superlight one, is then equivalent to 
requiring that the first eigenvalue Ai (or group of eigenvalues) be parametrically smaller than 
higher eigenvalues. It is interesting to note that there are general mathematical theorems 
which guarantee that such a hierarchy cannot occur if the compact metric 7 is Ricci-flat (or 
Ricci-positive). Indeed, if we consider, for simplicity, the Ricci-flat case, there are theorems 
(see [26], [27]) saying that there exist universal positive constants a n (d), b n (d) (which depend 
only on the dimension d of the compact manifold T) such that a n (d) 5~ 2 < X n < b n (d) 5~ 2 
for all n > 1, where 5 denotes the (metric) "diameter" of T. However, we wish to emphasize 
that, if one does not constraint the sign of the Ricci tensor, nothing prevents the occurrence 
of a spectral hierarchy. We conjecture that the generic situation where such a spectral 
hierarchy (between a finite group of abnormally small eigenvalues and the rest) occurs is 
a "near pinching" situation, i.e. the case where the manifold T is on the verge of getting 
split into two (or more) separate manifolds (of the same dimension d as T), as is illustrated 
in Fig. 3a. 

We have confirmed by some toy-model calculations that the near-pinching case (if the 
connecting "tube" between, say, two manifolds is not too long) does indeed lead to a 
spectral hierarchy. Let us also mention that a general theorem of Cheeger (see [28]) can 
be viewed as a (moral) confirmation of our conjecture. Indeed, this theorem says that a 
lower bound of the first eigenvalue is Ai > h 2 /A where Cheeger's constant h is defined as 
the lower bound of the ratio |S'|/inf (|Ti | , |r 2 |) when S runs over all closed submanifolds 
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a) b) 

Figure 3: Manifold T (Fig. 3a) is on the verge of splitting into two classically disconnected 

manifolds Ti and T 2 (Fig. 3b). These two manifolds may be connected at the quantum 
level. 



of T (of dimension d — 1) which partition V into two open manifolds Ti, r 2 , with common 
boundary dTi = dT 2 = S. We use the notation |T| to denote the (riemannian) volume of 
T. Note that h has the dimension of an inverse length, and that "pinching" does indeed 
correspond to the case where h — > 0. 

Physically, we can view the very light mode arising in a nearly pinched configuration 
r x U {tube} U T 2 as coming from the effect of a weak coupling between two "resonators" 
(or quantum mechanical systems) having regular spectra A^, A^ 2 ). Before coupling, the 

ground state is degenerate, A^ = A^ = 0. Weak coupling is generally expected to split 
this degeneracy into a doublet. As A = is always an exact eigenvalue of the combined 
system (corresponding to (po(y) = 1), this mechanism always leads to a small Ai (going to 
zero with the coupling). Note that the eigenmode corresponding to Ai is approximately 
equal to <p\{y) ~ e(y) where e(y) is +1 over T\, and —1 over r 2 . 

We are aware of the fact that weakly coupled string theory suggests compactification 
on Ricci-flat manifolds (which exclude a spectral hierarchy) . However, we think that string 
theory might still, in certain circumstances, allow for a spectral hierarchy: either because of 
a'-corrections to Einstein equations (which lead one away from the Ricci-flat case), or (more 
speculatively) because of conceivable quantum tunnelling effect between two (separate, but 
"near") Ricci-flat manifolds T\, r 2 . Pictorially, such a tunnelling situation is the limit of 
Fig. 3b where the link between Ti and T 2 is classically broken. The exponentially small 
coupling associated to such a tunnelling situation would naturally induce an exponentially 
small Ai, and thereby a bigravity coupling scale ji exponentially smaller than the string 
scale. 
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3.3 Bigravity and Connes' non-commutative geometry 

Within his general non-commutative geometry programme [29], Connes introduced the 
model of a two-sheeted space X, made from the product of a continuous space Y by 
a discrete "two-point space" {a,b} (or Z 2 ) : X = Y x Z 2 . Though the algebra A of 
"functions on X" (defined as the algebra of pairs of functions viewed as diagonal matrices 
diag(/ a (y), fb(y)) with y e Y) is commutative, the bimodule of 1-forms on such a space 
is not commutative [29]. Generalizations of this model (also based on the product of a 
continuum by a discrete space) were used in [30] to give a geometrical explanation of the 
structure of the Standard Model. In particular, it was found that the VEV of the Higgs 
field is related to the (non-commutative) "distance" between the two sheets. The metric 
aspect of such a two-sheeted space was developed along different lines by several authors 
[31], [32], [33], [34]. For instance, Ref. [31] introduced (non- commutative) analogues of the 
Riemannian metric, curvature tensor and scalar curvature, which enabled them to introduce 
a generalized Einstein-Hilbert action. This generalized Einstein-Hilbert action was found 
to contain (besides the standard integral of the scalar curvature of Y) a minimally coupled 
massless scalar field a related to the "distance" between the two sheets by d oc e~ a ( y \ 
An alternative approach to studying gravitational effects within general non-commutative 
spaces has been proposed in [34]. We shall follow this approach which is based on a general 
"spectral action principle" . In its simplest form, this principle is proposing to take as bare 
bosonic (Euclidean) action for any non-commutative model X the trace of the heat kernel 
associated with the square of the (non-commutative) Dirac operator D x of X: 

I = Tr exp (-*£>£■) • ( 50 ) 

Here t = m 2 introduces a cut-off, roughly equivalent to keeping only frequencies smaller 
than m . The cut-off-dependent Euclidean action (50) is viewed (a la Wilson) as the bare 
action at the mass scale ~ m . 

It seems that all previous works interested in the metric aspect of a two-sheeted space 
X = Y x Z 2 have restricted themselves (either for simplicity, or because of some constraints 
[33]) to the case where the metric is the same on the two sheets. By contrast, we focus 
here on the case where the two metrics are different, say g^ u and g^, and the aim of this 
subsection is to compute the "potential" V{gL,gn) implied by the spectral action (50). 
Following Connes (see p. 569 of the English edition of his book [29]) we define a Dirac 
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operator on a bi-Riemannian space X = Y x Z 2 as 



/ PL 



7 5 m \ 



D x = 



(51) 



\ 75 m 



Pr J 



This operator acts on bi-spinors (J},^) living on Y x Z 2 - Our conventions are that the 

(Euclidean) gamma matrices are hermitian, as well as 75 (which satisfies 75 = 1 and which 
anticommutes with the gamma matrices and therefore with the separate Dirac operators 
p L and Pr). The explicit form of the (hermitian) Dirac operators on each sheet is 



The explicit form of the spin connections f2 M will not be important for our calculations. 
On the other hand, the explicit form of the gamma matrices will be crucial. They read 
1l = 1° 1 EaLi 1r = l a EaR where {7 a ,7 6 } = 2 5 ab is a standard set of (space-independent) 
gamma matrices and where E% L (x), E% R (x) (where x G Y) are vierbeins corresponding 
to the two positive definite metrics g^ u , g^ u given on the abstract manifold Y : g^ u (x) = 
8 ab E^ L E^ L) etc. Note that the structure (51) assumes that we are given not only an 
identification map between corresponding points of the two sheets (here gauge-fixed by the 
identification of the two underlying abstract manifolds and the use of only one coordinate 
system x^ to describe the metrics on the two sheets), but also a one-to-one map between 
the spin structures, and in particular between any choice of vierbein. In other words to 
any E^ L must correspond a unique E^ R so that an arbitrary, local SO (4) rotation of E^ L 
corresponds to the same rotation of E% R . It is most natural to use as map E% L — > E% R the 
canonical map defined in [35] . This map can be defined by requiring that it reduces to simple 
rescalings E% R = e _Ma//2 E^ L when considering bi-orthogonal frames (as in Eq. (9) above). 
For simplicity, the quantity m in Eq. (51) (which "connects" the two sheets) will be taken 
to be a constant real scalar. More generally, it could also be a matrix when considering 
multiplets of fermions and could be space-dependent. We shall see that m is connected 
with the coupling scale fi in Eq. (6). It might be interesting to consider generalized models 
where m (and therefore jj) is linked to a fluctuating scalar m(x) oc e~ a ^ as in [31]. 



p L ,R = tit R (d, + n,(L,R)), Wni} = 2gr, i = l,r. 



(52) 
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The square of the Dirac operator (51) is easily obtained as 

/ pl + m 2 m l5 (p R -p L )\ 
\m l5 (p L -p R ) fll + m 2 ) 



D 2 - 
u x — 



(53) 



The heat kernel expansion of (50) is a series in increasing powers of t = m 2 which 
starts at order t~ 2 = m„. At this leading order the action I leads to two bare cosmological 
constant terms mQ/(47r 2 ) / d 4 x(^/g~Z+ \fg~R) ■ ^ ^ ne nex ^ t° leading order, 0(t~ r ) = 0(i71q), 
one gets two separate Einstein actions mg/(487r 2 ) / d A x(—^g~[ Rl — ^/cirRr) (with negative 
signs, as is appropriate for an Euclidean action Ie which is essentially Ie = —-f Minkowski with 
a positive signature metric) as well as a "potential" term + / d A xV(gL, g R ) proportional to 
m„m 2 . In view of its m^m 2 scaling, the potential V contains no derivatives of g L or g R . It 
can therefore be evaluated by considering the case of constant metrics g^ v , g R . We can then 

neglect the spin connections in (52) and go to the momentum representation (i <9 M — > k^) 
to set 



n 2 _ 
u x — 



/ fc 2 L + m 2 
V rwj^L -fa) 



(54) 



m-y 5 (Jc R -fa) \ 
/z 2 R + m 2 

Here jk L = 7^ ]k R = 7^ Using the explicit vierbein expressions of 7^ and 7^ we can 
rewrite these as fa = ^ a k^,fa = ^ a k R , where 



hi 



ft 



(55) 



In terms of these two different vectors (that live in a local Euclidean space common to the 
tangent spaces of the two sheets), one easily finds that the eigenvalues of D 2 X are 



An 



ki 



+ m 2 ± 



ffcl 



k R ) 2 



+ m 2 (k L - k R ) 



(56) 



Here, all squares are evaluated with the flat Euclidean metric 8 ab appropriate to the local 
Euclidean space where both k% and k^ live. In the limit m 2 — > (appropriate to the heat 
kernel expansion) the eigenvalues (56) read (we henceforth suppress the boldfacing of the 
Euclidean vectors k L and k R ) 



A, 



1 2 1 2 1 2 
k L + m + m 



(kL - k R ) 2 
h 2 — h 2 

n>T. ">\ 



+ 0{m 4 
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A_ = k\ + m 2 - m 2 l , L 2 7/ + G(m 4 ) . (57) 



h 2 — h 2 



The heat kernel action reads 



/ = Tr exp {-tD 2 x ) = 4 J >x / ' -J^ [e"< A + + e~ tx -}, (58) 

where the 4 comes from the trace in spinor space and where d i k is the fourfold integral over 
the covariant components k^. Expanding (58) in powers of m 2 leads to the mixing term 
+ 4m^m 2 / d A x (Vi + V 2 ) where 

^^/SW^-^' (59) 



V, = -t 2 j^-{e-^ + e~^). (60) 



d 4 A: 
(2^ 

After our factorization of wig = in front of V\ and V 2 , the expressions (59), (60) are 
easily seen to be t-independent. We can then evaluate them by setting, say, t — 1 in them. 
Noting that fc| = <5 afe A;^ k^ = g£ u k^ k v , etc., V 2 is easily evaluated: 

V 2 = -c 4 ((det + (det^)"3) = -c 4 ( v ^I+ (61) 



where c 4 = (levr 2 )^ 1 and where g L = det g^ u = (det^") -1 . The potential Vi is much more 
tricky. However, it can be nicely expressed by introducing a Schwinger-type parameter a 
(varying between and 1) and by using the identity (e _a — e~ 6 )/(a — b) = — Jq 1 daexpf— ((1 — 
a)a + a&)]. This naturally leads to the introduction of a one-parameter family of metrics 
g{a) interpolating between gi and gn (reached, respectively, when a = and a = 1). More 
precisely we define 

g^(a) = (l-a)gff + ag^. (62) 

[Note that the "line" connecting g^ to gn is "straight" when expressed in terms of con- 
travariant metrics (which naturally appear in the squared Dirac operator D 2 = —g^ = 
+ g^ v k^kv), but will become "curved" when expressed in terms of the covariant compo- 
nents g^ v (a)] In terms of the definition (62) (and the associated g^ipt) g av {pi) = 5^, 
g(a) = det g fll/ (a) = (det ^"(a)) -1 ) we find that Vi can be written as 

V, = + l - c 4 AE£ AE: £ da yfgjajg„ u (a) , (63) 
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where AE% = E^ R - E% L . Finally, the 0(m 2 ) piece of the Euclidean action (i.e. the "po- 
tential", remembering that £e uc i. = -^Mmk.) predicted by the non-commutative approach 
to two-sheeted spaces reads + / d 4 x V with 



= m 2 m 2 
' ' 4vr 2 



X - AE% AE v a jf* da Jgjaj g^(a) - ^g~l - y/gR 



(64) 



The explicit evaluation of the a-integral in (64) can be reduced to (incomplete) elliptic 
integrals. In fact, it can be reduced to the evaluation of the single integral 



I(9l,9r)= f daJg(a) = [ - d ° = (65) 

by using the identity 

( d d \ If 1 I 

I + J 7 (^' 9r) = ~2J da y9( a )9^(a). (66) 

\ 9 L 9 R / ^ 



When considering a bi-orthogonal frame, say with E^ L = diag(e ^^ 2 ), E^ R = diag(e ^/ 2 ) 
(so that g^ v = diag(e+^), g* v = diag(e+' 4 «) and g^ 1 g R = diag(e^) with fi a = n*-n%) the 
integral (65) is a rather simple elliptic integral of the first kind which, in principle, can be 

explicitly expressed in terms of the eigenvalues \^ ,R = e fla ' R . Of more direct interest for us 
is the discussion of the "weak-excitation" limit of V, i.e. the limit \x a = /x^ — fi^ — > 0, i.e. 
(9I 1 9r)u = K + h% with h% — > 0. In this limit we find that V behaves as (with a± — E /i a , 
a 2 = E /x 2 as above) 



2 2 
wig m 

4vr 2 



)'(- 2+ f-i) 

(g L g R )i(-2 + lhlK-^(K)*) . (67) 



Besides a negative m 2 -dependent, contribution to the cosmo logical constant (which has 
anyway bare contributions 0{mfy), we see that we do not get a Pauli-Fierz-type mass term 
for weak excitations away from = g R . We get instead (remembering that the (bare) 

Planck mass is M| = M R = ml/ (487r 2 )) a mass term proportional to m 2 h^h 1 ^ — \ (h 1 //) 2 
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Such a mass term contains a scalar ghost, but has the virtue (contrary to the Pauli-Fierz one 
~ m 2 [h u ^ h^ — (h^) 2 }) of exhibiting excellent continuity properties of the limit m 2 — > for all 
processes linked to the generation of gravitational fields by sources (see, e.g., Appendix C 
of [37] where it is easily seen that a — 1/2 leads to an Einstein- like propagator h^ v T^ u oc 

T^D-m 2 )- 1 (T^-\^ U T)). 

4 Phenomenology of Multigravity 

4.1 Bigravity, and bicosmology, versus massive gravity 

There is quite a sizable (and somewhat confusing) literature about the "problems" raised 
by having either "massive gravity" (i.e. a kind of finite-range version of Einstein's theory), 
or a "massive graviton" in addition to Einstein's massless one. We leave to a future publi- 
cation a detailed discussion of such issues, but wish to emphasize the fact that the change 
of paradigm, brought by focusing on a fully nonlinear bigravity theory, drastically modifies, 
in our opinion, the way one should view the traditional "problems" of massive gravity (in 
both senses recalled in the sentence above). One of the basic points is that many of the 
"problematic" issues (such as, unboundedness of the energy, singularity of the infinite-range 
limit) simply loose their meaning in a general bigravity setting. Indeed, these problematic 
issues make sense only for states (in some theories) which are, at least asymptotically, 
close to some trivial, Poincare invariant background. We think that, even when consider- 
ing formally "small" excitations above a trivial background state g^ u (x) = g^ u {x) = rj^, 
the exact bigravity configurations will generically develop into full-blown "bi-cosmological" 
configurations with fields that grow so much (in time and/or in space) so as to be outside 
the usually considered domain of bi-asymptotically flat configurations containing localized 
excitations. Note that most of the results concerning the "discontinuity" of the m 2 — > 
limit [38], [39], [37] implicitly (or explicitly) assumed such a framework of asymptotically 
decaying perturbations of a (minimum energy) Poincare invariant background. We think 
that, if one relaxes this asymptotic restriction, there exists a sector of bigravity theories 
which exhibits "physical continuity" for small m 2 , at the cost of cosmological behaviour 
on large scales. Note that such a claim, while being consistent with the works [40, 41, 42] 
which found continuity of massive graviton interactions in maximally symmetric ((A)dS) 
cosmological backgrounds, is somewhat different from the claim of [43], [44]. Indeed, the 
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latter claim seems to insist on a framework (and a language, like that of propagators, cou- 
pling and scattering states) which preassumes the restriction to localized excitations of a 
Poincare-invariant vacuum, i.e. that the metrics under consideration are asymptotically 
flat. Leaving to a future publication a detailed discussion of the "discontinuity" issue, we 
shall content ourselves here to sketch the general dynamical structure a la Arnowitt-Deser- 
Misner (ADM) [36] of bigravity theories. 

4.2 ADM analysis of bigravity theories 

We consider a general bigravity action (6). Let us decompose the two spacetime metrics 
4 4 #il into the two la P ses n l, N R (N L = (- A g L °)^), the two shift vectors N* L , N* R 
( 4 9ij n l = 4 0m) and the two spatial metrics gf p g§ (gf- = A g^ j ). We have 

ds 2 L = -N 2 L dt 2 + g^[dx l + N[ dt)(dx ] + N J L dt) , 

ds 2 R = -Nl dt 2 + gf^dx 1 + dt) (dx j + N j R dt) . (68) 

After integration by parts, each separate "Left" or "Right" pieces of the action (6) reads 
(say for the Left piece) 

I L = JdtJ d 3 x g% + II L $ L - N L H L - N{ H?) (69) 

where ir^ is the Left gravitational momentum density, Ul is a (generic) matter momentum 

density and where the left super-Hamiltonian, and super-momentum densities have the 
structure 

H L = ^^= {*t 4 -\*£)-Ml^g- L R L + Hr ttev , (70) 



M 2 L Jgi 

Hf = -2 Df ir j Li + H^ matter . (71) 

Let us now consider the interaction term — ( 4 (?l 4 fi , j?)^ V { a 9l 1 A 9r)- Using the fact that the 
local scalar V must be (in particular) invariant under transformations of the type dt' = A dt, 
dx' 1 = Aj {dxi + dt) one finds that it can only depend on the lapses and shifts through 



the combinations Nr/Nl and (N R — N l L ) / \J N L N R . Let us then replace the 8 variables N L , 
Nl, Nr, N r by the combinations 



/AT .1 ATI ATI 

N = \ N R N L , n = J— ,N l = -(N i R + Nl) , n* = R . (72) 

V H ' " Nl 2 v R L!l 2 VNrN l 
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With these definitions it is found that the total action reads 

1 = Jdtjd 3 x (tt£ g% + ix% g* + n L $ L + U R $ R - H) (73) 

where the total Hamiltonian density reads (here and below g^, gfj are the spatial metrics 
and g L = detg^, g R = detc/g) 

H(N, 7T, n, n\ g L , tt l , g R , n R , $ L , Tl L , $ R , U R ) ^H + jrH, (74) 

where 

H(n, n\ g, tt, $, n) = n~ x H L + nH R -n l + n l H? + (g L g R )^ V(n, n\ g L , g R ) , (75) 

H l (g,7r,$,Il) = HL + H t R . (76) 

The crucial point for the present discussion is the separation of the 8 lapse and shift 

variables into two sets: (i) the four "average" lapse and shifts N, N\ which are true 
Lagrange multipliers appearing only linearly in the action, and (ii) the four "relative" lapse 
and shifts n, n % which enter algebraically in the action (no kinetic terms) but in a non linear 
manner. The four average lapse and shifts give rise to four constraints, which are linked to 
the symmetry of the action under common diffeomorphisms: 

H = 0, Hi = 0. (77) 

N, N* are gauge variables which can be gauged away (e.g. to N — 1, Ni — 0). The four 
(first-class) constraints (77) can be used, together with the field equations (which involve iV 

and N l ) to eliminate four degrees of freedom (i.e. eight functions of positions and momenta). 
By contrast, the four relative lapse and shifts n, n l are not (undeterminable) gauge variables 
but are dynamical variables which are instantaneously determinable in terms of the other 
variables (g,ir, n) by their (algebraic) equations of motion: 

dH dH 

-- = 0, 7r- = 0. 78 

on on 1 

This result generalizes the findings of [37] which studied the case of "massive gravity" , i.e. 
(21). We must assume here that the potential V has a "good" dependence on n and n l 
which allow for an (essential) unique solution of Eqs. (78) for a generic (or at least an 
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open) domain of free dynamical data g, ir, <3>, n. We think that the only (covariant) sit- 
uation where n and rii combine with N, TP to generate more (gauge-related) Lagrange 
multipliers is the case where V is linear in Nl and Nr, which must then correspond 
(by covariance) to V = c L N L ^gJ + c R N R ^g^ = c L (-det 4 ^) 1/2 + c fl (-det 4 ^) 1/2 . 
For instance, we can think that V(n, rii) contains terms quadratic in ri l (as already fol- 
lows from a Pauli-Fierz mass term), and behaves, for both large (respectively, small) 
n as + an 2 (respectively +bn~ 2 ). Note also that if we define the new scalar poten- 
tial V ncw by factoring (-det 4 ^) 1/2 + (-det 4 £*) 1/2 instead of (W) 1/4 from V, 
i.e. V = ((-det 4 ^) 1/2 + (-det 4 gj? u ) 1/2 )V new , the last term in Eq. (75) will become 

{ n ~ 1 y/9~L + ^i/^)^™^' Ui ' 9 L > 9 R )- ^ * s then enough to require that V new (n) grows in 
any manner (even logarithmically) towards +00, as n — > +00 or n — > 0. Such conditions 
ensure the existence of (possibly non unique) solutions of the equations of motion of n and 
n\ Eq. (78). 

We can then use (78) to eliminate n and n % (by replacing them by their expression in 
terms of the other dynamical data). It is then easily seen that the reduced Hamiltonian 
H Te d(g, ft, n) obtained by inserting these expressions into (75) defines (together with 
(76)) a dynamical system for the variables g^, ir^, gf^ $ L , n L , $_r, n R (submitted 

to the four first-class constraints (77) coming from the Lagrange multipliers N, N*). For 
instance, if we consider the matter-free system, we end up with the 6 + 6 degrees of freedom 
linked to g^ and gf^ from which must be subtracted 4 degrees of freedom killed by the 
first-class constraints (77). This leaves us with 8 degrees of freedom. As in the analysis of 
(nonlinear) massive gravity in [37], which concluded to the presence of 6 degrees of freedom 
(instead of the expected 5 of a Pauli-Fierz linear graviton), we have here 2 + 6 = 8 (where 
the 2 can be formally thought of as corresponding to an Einstein (massless) graviton, and 
the 6 to a "massive graviton"). 

Two of the potential defects of the supplementary tensor degree of freedom (1 = 6 — 5) 
are, according to [37]: (i) the unboundedness of the total "energy", and (ii) experimental 
difficulties (e.g. with light scattering by the Sun), even if a suitable mass term can be 
found for which the m 2 — * limit exists. Our point of view concerning (i) is to argue 
that the notion of energy is not defined when considering (as we argue must be done) non- 
asymptotically flat metrics, with cosmological-type behaviour at infinity. Alternatively, we 
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can dismiss the problem of spatial boundary conditions by considering spatially compact 
manifolds (e.g. with toroidal topology). For such a situation, the dynamics associated to 
(73)-(78) should entail a well-defined (classical) evolution system for g L , ir L , g R , n R , . . . The 
ill-defined issue of "unbounded energy" is then transformed in a well-posed dynamical ques- 
tion: do Hamilton's equations of motion quickly lead to a catastrophic evolution towards 
some singular state?, or do they admit many solutions which evolve rather quietly on times 
scales comparable to the age of the universe? (which is the only stability property which 
is really required by experimental data). This question will be discussed in detail in [12]. 
Let us only mention here the result that there does exist, for suitable potentials V, many 
solutions which can quietly evolve on Hubble time scales or more. 

4.3 Phenomenology and a new form of dark energy 

Using the dynamical, and cosmological like, viewpoint expressed in the previous subsection, 
let us now briefly discuss why we think that bigravity is not only compatible with existing 
gravitational data, but might also furnish a natural explanation of the recently observed 
cosmic acceleration. Let us first argue that there exist large classes of bigravity data 
9Li^Li9Ri^Ri ■ ■ ■ which can adequately represent the universe as we see it at the present 
moment. For definiteness, we assume that we "live on the right brane" (when viewed in 
brane language), i.e. that the matter around as is made of <3>R-type matter only. Let us start 
by considering an instantaneous "Einstein" model of our universe, i.e. an exact solution 
gn, 7Tr, $r, Ur of the constraints H R = = Hf. Let us complete this configuration by 
a random "Einstein" model of the (shadow) left universe, i.e. a solution g L , ir L , $ L , Yi L of 
H L = = H[. Taken together, these two configurations "nearly" satisfy the bigravity 
constraints (77) and (78). More precisely, (77) is satisfied modulo a term proportional 
to V, while (78) is satisfied modulo terms proportional to dV/dn and dV/dn l . Let us 
assume that all dimensionless variables (n, ri 1 and gj} gn) are of order unity, and that 
V oc /i 4 is at most comparable to the average cosmological energy density V ~ 1CT 29 g cm" 3 
(i.e. n ~ 10 -3 eV) (in right units, say). Instead of viewing Eqs. (78) as equations for 
determining n and n l , we can pick rather arbitrary (initial) values of n and n l (or order 
unity) and slightly deform the Einstein data gn, n R , . . . , g L , til, ... to compensate for the 
small violation of the usual Einstein constraints brought by the terms proportional to V, 
dV/dn and dV/dn\ It is intuitively clear that there are many ways of doing so, i.e. of 
constructing exact bigravity initial data <?.r, 7Tr, 7Tl, .. . which exactly satisfy (77) 
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and (78) for arbitrarily given n and n\ Locally, say around our Galaxy, the new, deformed 
data g R , ir R , . . . can be constructed so as to be experimentally indistinguishable from a pure 
Einstein model (after all, we are simply modifying the stress-energy tensor in the Galaxy 
by less than lCT 29 gcm~ 3 , which is many orders of magnitude smaller than the average 
density in the Galaxy). If the dependence of V on n and n % is adequate the equations (78) 
will continue to admit a solution (n, n l ) during the future evolution of the other dynamical 
variables. In fact, as (under a general assumption made in Section 2 above) we know 
one exact (but physically trivial) solution of the full bigravity evolution equations, namely 
9i v it,x) = g$,(t,x), i.e. g L = g R , ix L = ir R , $ L = . . . with n = 1, n l = 0, we 
expect (by mathematical continuity) that there will be classes of bigravity solutions where, 
during a long time, g^ ~ g R , n L ~ ir R , . . . with n ~ 1, n l ~ 0. The crucial question is 
whether one can solve Eqs. (78) for a long time (without catastrophe) for more general 
data where gj} Qr = 0{1) and n = 0(1) = n\ This question will be addressed in [12] for 
cosmological-type solutions and in [13] for solar-system-type solutions. Note that this is 
here that the potential "discontinuity" problems linked to the m 2 — > (or /i — > 0) limit 
show up because the potential V(n,n l ) is proportional to /i 4 , so that, when solving for n 
and n l Eqs. (78), // will tend to appear in a denominator and might cause the solution n, n l 
to take parametrically large values, proportional to some negative power of ji (depending 
on the behaviour of, say, V{n) asn^+ooorn^O). 

Assuming, for the time being, the continuous existence of regular bigravity solutions, 
evolved from some data g^, ttl, g R , ir R , n, n l , . . . we can finish by mentioning some of the 
pleasing phenomenological aspects of bigravity. First, bigravity exactly satisfies the equiva- 
lence principle, because each type of matter (say whithin "our universe") is universally 
coupled to the corresponding metric, say g^ v . Second, (as just discussed) there are classes 
of bigravity solutions which differ from standard Einstein ones only by the presence "on 
the right-hand side" of Einstein's equations of numerically very small additional terms (say 
t% ~ ^.v ~ 10~ 29 gem 3 in the covariant form (14)), which locally modify g^ v and Q R (and 
9%ji om y m a numerically very small way (though they might globally forbid the stable 
existence of asymptotically flat models). These solutions will be fully compatible with all 
local (or quasi-local) experimental tests of relativistic gravity: such as solar-system tests 
and binary-pulsar tests. Third, if /i indeed happens to be of the order of 10~ 3 eV, and 
if g^ 1 g R is of order unity, bigravity will only lead to experimentally significant deviations 
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from Einstein's gravity on cosmological scales. Moreover, if, seen from our universe g^ v , 
we view g^ v as an "external field" , or, more precisely, if we (approximately) view the "dif- 
ference" between the two metrics gj} gn as a given (time varying) tensor "condensate" of 
order unity, the potential term V = ^JgHig^ gh)^ V^Ql 1 9r) can be approximately viewed 
as a time- varying "vacuum energy" term (of order //), i.e. as a kind of "dark energy". It 
is tempting to assume that this new form of dark energy (which might be called "tensor 
quintessence") can explain the observed cosmic acceleration. It might also be used in pri- 
mordial cosmological scenarios, possibly when using the idea mentioned above that \i could 
be an evolving field. See [12] for a study of this new form of dark energy, and its phe- 
nomenological differences with quintessence models based on evolving scalar (rather than 
tensor) condensates. 

5 Conclusions 

In this paper we suggested a new paradigm concerning "massive gravity" and "large scale 
modification of gravity" . Considering the fully nonlinear bigravity action suggests to change 
viewpoint: instead of the theory with massless and massive graviton(s) we had in linearized 
approximation, we are dealing with several interacting metrics. We introduced the concept 
of universality class which we formulated using bigravity (two interacting metrics) as an 
example. Different approaches (brane, KK, non-commutative geometry) naturally lead to 
different universality classes for the fully nonlinear bigravity action. Another important 
new suggestion is that almost all solutions must now be of the non-asymptotically flat 
(cosmological) type. 

This new formulation can change the standard problematic of the m 2 — > discontinuity. 
We showed the existence of classes of solutions that are compatible with "our universe". 
However, we do not claim to have proven that general solutions of bigravity are phenomeno- 
logically acceptable. The two main problems of massive gravity (ghost, potential blow up 
of some field variables when m 2 — > 0) must still be examined in detail. The important 
problem is to find the matching to the local sources of the field so that the full metric is 
free of singularities. We do not worry about matching at infinity because we abandon the 
requirement of asymptotic flatness. It is possible that in some models of bigravity such 
local matching does not exist because of the explicit or implicit presence of ghost modes in 
the theory. Such models would be physically unacceptable. We note in this respect that 
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the 6-dimensional model discussed in [21] which does not contain negative tension branes, 
contains instead either branes with equations of state violating the weak energy condition 
T brane ^ p, > Q ( with light _ like £m) or has a conifold singularity in the bulk. The physical 

consistency of this model must be further investigated. We have also quoted mathematical 
theorems linking the existence of a hierarchical spectrum (necessary for the derivation of 
an effective bigravity Lagrangian) to the necessary negativity of the Ricci curvature of the 
compactified manifold. This sign condition might hide the presence of ghost-like fields in 
the theory. These questions are pressing and deserve detailed investigation. 

Assuming a positive resolution of these issues or simply taking the phenomenological 
viewpoint that nonlinear bigravity Lagrangians open an interesting new arena for non 
standard gravitational effects, we shall explore in future publications [12], [13] the nonlinear 
physics of bigravity actions, with a particular view on its cosmological aspects, as it may 
provide a natural candidate for some new type of "dark energy" . 
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Appendix 

In this Appendix we check the consistency of the linearized limit of the nonlinear action 
(6) with a direct linearized analysis of the coupling strengths of massless and light graviton 
modes in brane models. Omitting the tensor structure (and considering only the relative co- 
efficients between the various terms) the Lagrangian describing the coupling of the massless 
graviton mode h°, and of the lightest one h 1 , reads 

L Lm = h°d 2 h° + h 1 (d 2 + mj) h 1 + ^ f (h° + aLh 1 ) T L + ^- (h° - anh 1 ) T R (79) 

IVlp IVlp 

where the coefficients oll,r describe the relative strengths of the massive graviton coupling 
to the matter on left and right branes. It seems that there are four parameters here: 
M p , mi, ocl, and or. But actually = a^ 1 which is extremely important as we shall 
see next. This relation follows from the expression for a R which was obtained in [1] (see 
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Figure 4: Here L\ is the position of the bounce. The left configuration is just the mirror 
image of the right one and the positions of the bounces are related by L\ = L 2 — L 1 . Under 
this transformation left and right branes are exchange their roles and, at the same time, 
a R = exp(2Li - L 2 ) —> a R = exp(2L 2 - Li) = a^ 1 . 



Eq.(20) and Eq.(22) there) 

a R = exp (2Li - L 2 ) . (80) 

Being derived originally for the H h model this expression holds for other models with 

bigravity, for example the ++ model. In Figure 3 it is shown that one can interchange left 
and right branes by changing the position of the bounce from L\ to L\ = L 2 — L\. One 
gets the new coupling strength 

a R = exp (2Li - L 2 ) = exp (L 2 - 2L X ) = — . (81) 
v ' ol r 

At the same time it is easy to see that a new right brane is just an old left one, so that we 
have the result 

oil = olr = — . (82) 
ot R 

This relation is crucial to the consistency of the nonlinear bigravity approach because only 
in this case can one relate ho and hi to hi, ~ (h° + ol/ii) and h R ~ (h° — a R hi) by 
orthogonal rotation. If it were not the case one would get mixing between hi and h R even 
in the limit mi = and we could not have two non-interacting worlds. Introducing 

hi = cosOh + sin Oh 1 , , h R = sin 0h° — cos Oh 1 , tan# = a R (83) 

we can rewrite (79) as 

L Lin = h L d 2 h L + h R d 2 h R + . /f2 mi . /r2 (M L h R - M R h L f + ^hiTi + ^rTr (84) 

M L + M R Mi M R 
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where M L = M p /y/l + a 2 , M R = aM p /y/l + a 2 

Let us note that in the limit L 2 — > 00 both ol and M p are divergent and — > 0. In 
this limit Ml is finite and Mr — > 00. The massless graviton /i^ becomes essentially a free 
sterile particle and decouples from the spectrum, while the massive graviton Hr interacts 
with matter on the left brane only. Long range gravity completely decouples from the right 
brane. 

Here we discussed the linearized bigravity lagrangian for flat branes, but one can get 
the same picture for (A)dS branes. The limiting case L 2 — > 00 with corresponds to a single 
AdSi brane was considered in [20] (see also [14, 45, 46, 47]). 
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